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ABSTRACT 

Possible detection of signatures of structure formation at the end of the 'dark age' epoch 
(z ~ 40 — 20) is examined. We discuss the spectral-spatial fluctuations in the CMBR temper- 
ature produced by elastic resonant scattering of CMBR photons on HD molecules located in 
protostructures moving with peculiar velocity. Detailed chemical kinematic evolution of HD 
molecules in the expanding homogeneous medium is calculated. Then, the HD abundances 
are linked to protostructures at their maximum expansion, whose properties are estimated 
by using the top-hat spherical approach and the ACDM cosmology. We find that the optical 
depths in the HD three lowest pure rotational lines for high-peak protohaloes at their max- 
imum expansion are much higher than those in LiH molecule. The corresponding spectral- 
spatial fluctuation amplitudes however are probably too weak as to be detected by current 
and forthcoming millimeter-telescope facilities. We extend our estimates of spectral-spatial 
fluctuations to gas clouds inside collapsed CDM haloes by using results from a crude model 
of HD production in these clouds. The fluctuations for the highest-peak CDM haloes at red- 
shifts ^20 — 30 could be detected in the future. Observations will be important to test model 
predictions of early structure formation in the universe. 

Key words: cosmology: first stars — galaxies: formation — molecular processes — cosmol- 
ogy: theory — dark matter 



1 INTRODUCTION 

In the last years a great deal of interest arose for understanding the 
formation of the first structures in the universe at the end of the 
so-called 'dark age' epoch. This interest is motivated not only by 
the possibilities of direct measurements of the physical conditions 
prevailing in these (proto)structures and the constraining of cosmo- 
logical parameters at very high redshifts, but also by the measure- 
ment of the primordial abundances of key elements (e.g., D and Li) 
in pregalactic epochs as a direct signature of the Big Bang Nucle- 
osynthesis (BBN). From the observational side, the re-emission or 
absorption of Cosmic Microwave Background Radiation (CMBR) 
photons in resonant lines from H2, LiH, and HD, among other pri- 
mordial molecules, have been suggested as a viable way to detect 
cosmic protostructures (Dubrovich 1977,1983; Maoli, Melchiorri 
& Tosti 1994; Maoli et al.1996). Resonant absorption by neutral H 
at its 21 cm transition has also been proposed to trace protostruc- 
tures at early cosmic times (Hogan & Rees 1979; for more recent 
studies see Barkana & Loeb 2005 and more references therein). 

The first attempt to detect LiH emission and Doppler-induced 
anisotropies in the CMBR from the protostructures at redshift 
z = 100 was made with the 30 meter IRAM radio-telescope (de 
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Bernardis et al. 1993). Only upper limits for the abundance of LiH 
were obtained. It was believed that the relative abundance of LiH, 
[LiH/H], at z = 30 - 100 could be rather adequate (up to 10" 10 ) 
to produce detectable spectral features due to the optical depth in 
the lines of the LiH rotational structure. However, more detailed 
calculations of the primordial LiH relative abundance have shown 
that, at the discussed epochs (z = 30 — 100), it should have been 
actually very small, [LiH/H] < 10~ 18 (Stancil, Lepp & Dalgarno 
1996; Bougleux & Galli 1997). Therefore the optical depth in the 
rotational lines of LiH in the primordial gas is too small to produce 
observational features (Galli & Palla 1998; Puy & Signore 2001; 
Maoli et al. 2004). 

The other primordial molecule, potentially interesting to de- 
tect features of precollapse structures from the dark epoch, is 
deuterated hydrogen, HD. The dipole moment of the HD molecule 
is actually much smaller than the one of LiH: cZhd = 8.3 x 10 -4 
debyes (Abgrall et al. 1982) and dLiH = 5.88 debyes. However, for 
[LiH/H] « 10" 18 and [HD/H]=(1 - 5) 10" 9 , one obtains in a first 
approximation thd/tlih « [HD]d^ D / [LiH] d^n ~ 20 - 100, 
concluding that HD is a species more suitable for detecting the pri- 
mordial structure signatures than LiH. Besides, the HD molecule 
is also able to amplify its emission due to ro-vibrational lumi- 
nescence effects (see e.g., Dubrovich 1997). At the moment, the 
most remote of the detected absorbers in rotational lines of HD 
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molecule is at redshift z — 2.3377 in the spectrum of the quasar 
PKS1232 + 082 (Varshalovich et al. 2001). Nevertheless, the 
emission of this molecule should be observable even from the pre- 
galactic epoch. 

In the intermediate epochs at 6 < z < 30, where nonlinear 
collapse takes place for the mass scales of interest, the HD molecule 
can reemit the kinetic energy of gas out of equilibrium. The emis- 
sion from primordial molecular clouds in the pure rotational struc- 
ture of HD molecule owing to cooling processes has been calcu- 
lated first time by Shchekinov (1986) and more accurate by Ka- 
maya & Silk (2003). Recently, Lipovka, Nunez-Lopez & Avila- 
Reese (2005) have revisited the cooling function of this molecule 
by taking into account its ro-vibrational structure. At temperatures 
above ~ 3000K and at high densities, the cooling function is 1 — 
2 orders of magnitude higher than reported before e.g., (Flower 
et al. 2000). If photo-ionization or shocks heat the dense collaps- 
ing clouds to temperatures above ~ 1000K, these new calculations 
show that the gas will be able to cool efficiently again with further 
HD resonant line emission, and the intensity of the collapse en- 
ergy reemission by the HD molecule calculated by Kamaya & Silk 
(2003) will increase probably due to the ro-vibrational structure of 
this molecule. Recently, it was also suggested that the H2 and HD 
fractions can be substantially large in relic HII regions (Johnson & 
Bromm 2006; Yoshida 2006), as well as in large mass haloes in 
which the gas is collisionally ionized (Oh & Haiman 2003). 

At the linear stage of evolution of the first structures of interest 
(z > 20 — 30), the elastic resonant scattering of the CMBR by pri- 
mordial molecules in the evolving protostructures can damp CMBR 
primary anisotropies or, if the protostructures have peculiar ve- 
locities, secondary anisotropies in emission (spectral-spatial fluc- 
tuations) can be produced (Dubrovich 1977,1993,1997). In these 
cases, internal energy sources are not required to produce the ob- 
servational spectral and spatial features. Therefore, primordial res- 
onant lines are ideal to study the evolution of cosmic structures 
during their linear phase or the first phases of their collapse, be- 
fore star formation triggers. The main goal of the present work is 
to study the possible observational signatures from these phases 
previous to the collapse of luminous objects in the Universe. 

The HD molecule have at least two advantages for this task: 
1) It is formed rather fast in the homogeneous primordial gas at 
z m 150, through channels similar to those of molecular hydrogen, 
and it may reach relatively high abundances at those epochs. In- 
deed, at z = 100, the relative abundance of HD is already ~ 10% 
of the presend-day one. 2) The redshifted wavelength range of the 
rotational structure transitions of the HD molecule is very comfort- 
able for detection. The first, second, and third rotational transitions 
at the wavelength of 112.1 /j,m, 56.2 p,m, and 37.7/otj can be de- 
tected at 0.8-3.5 mm if the lines are produced at redshifts around 
20-30, typical of the epochs when the first stars probably started 
to form according to numerical and semi-analytical results (for a 
review see Bromm & Larson 2004). Several facilities planned or 
under construction will cover the submillimetric and millimetric 
wavelength ranges with high sensitivity, for example the Large Mil- 
limetric Telescope (LMT, or GTM for its acronym in Spanish) in 
Mexico, whose radiometers cover the range from 0.8 to 3 mm. 

Galli & Palla (2002) presented results from detailed calcula- 
tions of chemical kinetics in thermally evolving primordial clouds, 
though they did not follow the hydrodynamical evolution. The rel- 
ative abundance of the HD molecule that these authors found is 
very high in the cooled gas clouds and, as it will be shown later 
here, must be observable, so, more accurate calculations of the 



protocloud and cloud formation process and emission by the HD 
molecule are required. 

The aim of this paper is twofold: to recall the importance of 
the HD molecule for investigating the very high-redshift universe, 
and to calculate the parameters of HD lines of interest. The latter is 
a crucial step for planning future observational strategies. We will 
estimate the optical depth for protoclouds in the rotational lines 
of the primordial HD molecule based on the A Cold Dark Matter 
(ACDM) scenario of cosmic structure formation. For alternative 
scenarios (e.g., Khlopov & Rubin 2004), the predictions will prob- 
ably be different. 

With the aim to obtain useful estimates, here we calculate de- 
tailed chemical kinematics evolution for HD molecule in the ex- 
panding homogeneous medium (§2). Then we calculate the prop- 
erties of high-density ACDM mass perturbations (protohaloes) at 
their maximum expansion (§3.1), and estimate the optical depths 
in the three lowest rotational line transitions of the HD molecule 
fraction in these protohaloes (§3.2). We bear in mind that these 
opacities will slightly decrease in the first phases of collapse when 
the line-widths increase and become potentially observable as sec- 
ondary spatial - frequency anisotropies in the CMBR. We estimate 
the amplitudes and angular sizes of these anisotropies assuming the 
linear model for the peculiar velocities of the mass perturbations, 
and explore the possibility to detect them by the LMT/GTM sub- 
millimeter telescope under construction in Mexico (§3.3). Finally, 
we present a discussion of the implications of our results (§4). 



2 EVOLUTION OF THE HD MOLECULE ABUNDANCE 

Following, we calculate the evolution of the HD abundance in the 
homogeneous pregalactic medium. This is a good approximation 
to the HD abundance in mass perturbations at their linear evolu- 
tion phases as yet. We use the concordance cosmological model 
(ACDM) with fi M = 0.04, « M ,o = 0.27, ft A = 0.73 and 
Ho = 71 kms" 1 Mpc" 1 (e.g., Spergel et al. 2003). Standard BBN 
yields for the light elements are assumed. 

In the early epochs, before the first cosmic objects are formed, 
the thermodynamics of the primordial gas is out of equilibrium. 
Therefore, in order to describe the molecular dynamics, we need to 
solve the kinetic equation system, which can be written as follows: 

—j^- = njrikRjki — Tli n m Rimn + 

j,k m.n 

^^UjRji — Ui'^^Rim, (1) 

3 tn 

where t is the cosmic time, m is the density of the species i, 
Rjki{Tk) are the rates of the collisional processes j + k — > i , 
as functions of the kinetic temperature Tk, Rji(T r ) are the rates of 
the radiative processes (formation and destruction of the molecule 
by the CMBR photons characterized by the radiative temperature 

T r ). 

In eq. Q the positive and negative terms correspond to the 
formation and destruction processes of the species rii, respectively. 
To solve the equations system Q, it is convenient to introduce new 
variables: Xi = m/nB, where tib = ns,o (1 + z) 3 and ns,o is 
the present-day baryon number density, ub.o = fisn cr . It is also 
convenient to substitute the cosmic time variable with the redshift, 
by using the appropriate relation for the flat cosmology with cos- 
mological constant adopted here: 
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dt = 



dz 



H o (1 + z) 0WlW+ftA 
With these new variables, eqs. Q take the form: 

dxi 1 



(2) 



dz Ho(l+z) ^/fiM,o(l + z) 3 + A ' 

«0 (1 + «) 3 I XjXkRjki — Xi ^ XmRimn 



(3) 



As was mentioned above, the reaction rates Rijk and i?^ are 
functions of the radiative temperature T r and the kinetic temper- 
ature Tk, which depend on many parameters. The radiative tem- 
perature can be taken as in the case of the adiabatic expansion, 
T r — To (1 + z), but in order to find the kinetic temperature we 
have to solve the differential equation for Tk : 

dt 



^ = -2T k H a ^/n M , (i + zr + n A + 



+ 



8a t abT r 



X e (T r — Tk); 



(4) 



3m e c 

where the first term in the right-hand part describes the change 
on the temperature due to the universe expansion, the second one 
corresponds to the thermodynamics of the primordial molecules, 
and the third one is Thomson scattering. The second term can be 
written as follows: 

'dtk\ 2 (r mo ; — A mo j) ^ Q ch ^ 2T k /dn\ 

\dt J ch 



/dtk\ 

V dt I m ol 



(5) 



oi ink 3nk 

where Y mo i and A mo i are the heating and cooling functions, <d c h is 
the energy gain (loss) due to the chemical reactions, characterized 
by their enthalpy AH, 



reactions 



RijknirijAH, 



(6) 



The summation is carried out over all chemical reactions. 

The factor ) in the last term of eq. corresponds to the 
change in the primordial gas density due to the chemical reactions, 
which do not conserve the number of initial species (for example 
A + B^C, or A + B ^ C + D + E). 

The calculations of the molecular abundances were carried out 
for a wide range of redshifts for the following species involved in 
HD molecule formation: H, H+, Hj, H 2 , H - , e", H+, He, He+, 
HeH+, D, HD, H 2 D+, D+, D , and HD+. Other species are neg- 
ligible in the formation processes of the HD molecule due to their 
small abundances. The following reactions are of particular impor- 
tance in the formation/destruction processes of the HD molecule in 
the primordial gas: 

D + H 2 = HD + H 
D+ + H 2 = HD + H+ 
D~ + H = HD + e" 
H~ + D = HD + e" 
H + HD = H 2 + D 

H+ + HD = H 2 + D+ (7) 

More detailed consideration can be found in the many papers (see 
for example Galli & Palla 1998,2002). The results of our molec- 
ular abundance calculations are shown in Fig. 1. One can see that 
the relative abundance ratio [HD]/[H 2 ] at low redshifts is slightly 
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Figure 1. Fractional abundance of chemical species involved in the HD 
molecule formation as a function of redshift for the expanding and adiabat- 
ically cooling homogeneous primordial medium. The concordance ACDM 
cosmological model was used. 



larger than 10 , whereas at high redshifts this value tends to the 
assumed primordial ratio [D]/[H]=3 10" °. This fact is explained 
as the deuteration process, D + +H 2 — > H + +HD and D + H 2 — » 
HD + H, that leads to relatively high abundance of the molecule, 
up to [HD]/[H 2 ]= 10" 2 ' 39 at * = 10 in collapsed clouds (see e.g., 
Galli & Palla 2002). 



3 THE HD OPTICAL DEPTH IN PROTOHALOES AND 
HALOES AT HIGH REDSHIFTS 

We are now interested in obtaining estimates related to the observ- 
ability of the line features of the primordial HD molecule in proto- 
structures in their linear or quasi-linear evolution regime, i.e. be- 
fore gravitational collapse, and first star formation. As mentioned 
in the Introduction, previous works have shown that the largest res- 
onant line opacities in homogeneous mass perturbations are pro- 
duced when they are at their maximum expansion or turn-around 
epoch. Based on these conclusions, here we will calculate the HD 
resonant line opacities at the turnaround epoch of CDM high-peak 
mass perturbations. After the first bursts of star formation (Popula- 
tion III stars) in the high-peak (3 — 6a) collapsed perturbations, re- 
ionization, feedback and metal enrichment processes change drasti- 
cally the properties of the medium in the lower density mass pertur- 
bations. This is why we will focus on the protostructures emerging 
from high 3-6<r peaks. The interaction of CMBR photons with the 
primordial molecules inside these protostructures will produce sig- 
natures related to the pristine conditions of gas before the first star 
formation in the universe. 



3.1 Properties of 3 and 6a top-hat CDM protohaloes and 
haloes 

To calculate the optical depth of primordial protoclouds by scat- 
tering of CMBR photons on HD molecules we need a model for 
these protoclouds. We use the spherical top-hat approach to calcu- 
late the epoch, density, and size of CDM mass perturbations at their 



4 Nunez-Lopez, Lipovka & Avila-Reese 

Table 1. Properties of ACDM mass perturbations 
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maximum expansion for the ACDM cosmology adopted here. This 
approach is justified for the kind of estimates in which we are inter- 
ested (e.g., Tegmark et al. 1997), and for the quasi-lineal phases of 
gravitational evolution. It is assumed that gas is equally distributed 
as CDM in the mass perturbation but with a density £Ib,o/Qm,o 
times lower. According to the spherical top-hat approach, the over- 
density of a given mass perturbation, 5(M), grows with z propor- 
tional to the so-called growing factor, D{z), until it reaches a (lin- 
early extrapolated) critical value, S c , after which the perturbation is 
supposed to collapse and virialize at redshift 2 co i (for example see 
Padmanabhan 1993): 



5(M;z col ) = 8 (M)D(z co i) = 6 B ,o- 



(8) 



The convention is to fix all the quantities to their linearly extrapo- 
lated values at the present epoch (indicated by the subscript "0") in 
such a way that D(z — 0) = Do = 1. 

The redshift of turn-around, z ta , is calcualted by using 
the top-hat sphere result of ita = 0.5t co i • Following the evo- 
lution equation for the top-hat sphere, one finds for our flat 
ACDM cosmology that the average sphere density at zt& is 
given by p ta ~ 5.6pb g (zta), i.e. similar to the Einstein-de Sit- 
ter case (Wang & Steinhardt 1998; Horellou & Berge 2005). 
The background density at z ta is Pbg(zta) = 1-88 10 -29 g 
cm -3 n_M,o'i 2 (l + Zta) 3 - The corresponding numerical gas density 
is n t a = {p%a,/m p )(Q.Bfi/^M,o), where m p is the H atom mass. 
Finally, the radius of the sphere of mass M at the turn-around or 
maximum expansion is given by 



Rta, 



3M 

47rpta 



(9) 



We need now to connect the top-hat sphere results to a char- 
acteristic mass of the perturbation, M, produced within the ACDM 
structure formation scenario. The primordial fluctuation field, fixed 
at the present epoch, is characterized by the mass variance, a(M), 
which is the rms mass perturbation smoothed on a given scale R 
corresponding to the mass M. For our calculations, we use the 
power spectrum shape given by Bardeen et al. (1986) normalized 
to erg = 0.9. Thus, the mass perturbation overdensity linearly ex- 
trapolated to z — can be defined as: 



where M and v are the perturbation mass and peak height, respec- 
tively. For average perturbations, v = 1, while for rare, high- 
density perturbations, from which emerged the first structures, 
v » 1. 

By introducing eq. 1 101 into eq. {8j one may infer z co u and 
therefore z t a, Pta, nta., and R ta , for vo mass perturbations from 
the ACDM primordial fluctuation field. In Table 1 we present these 
quantities for 3<j perturbations of different masses. As mentioned 
above, we are interested in HD line signatures from the protoclouds 
before first luminous objects in the universe formed and started to 
re-ionize it, and therefore only high-peak low mass perturbations 
will be considered. 

Cosmological N-body + hydrodynamics simulations showed 
that the first stars may have formed indeed at redshifts as high as 
~ 18 - 30 in (~ 3 - 4<r) CDM mini-haloes of M~ 10 6 M Q (Abell 
et al. 1998; 2002; Fuller & Couchman 2000; Bromm et al. 2002; 
Yoshida et al. 2003). Related to this is the fact that haloes of any 
given mass that collapse first (high peaks) do not populate "typical" 
regions at all, but rather "protocluster" regions (White & Springel 
2000; Barkana & Loeb 2002). In a recent paper, Gao et al. (2005) 
used a novel technique to follow with unprecedentedly high reso- 
lution the growth of the most massive progenitor of a supercluster 
region from z ~ 80 to z — 0. By using this simulation, Reed et 
al. (2005) have found that, when the mass of the progenitor halo 
is w 2.4 lO 5 h -1 M0 at z w 47, it undergoes baryonic collapse via 
H2 cooling, triggering star formation at this early epoch. Of course, 
this halo should emerge from a very rare fluctuation peak, v — 6.5 
(Reed et al. 2005). Notice that rare peaks are strongly clustered, so 
that the probability to find a cluster of these (proto)halos is high. In 
the second part of Table 1 we present the same as in the first part, 
but for 6(J haloes. 



3.2 Optical depths 

The optical depth for a protocloud of size L is given by: 



tv(L) 



tt v (x) dx, 



(11) 



<So(M) = fff(M), 



(10) 



where a„ (a;) is the absorption coefficient, and integration is car- 
ried out over the line of sight. In the case of a mass perturbation at 
its turn-around epoch, when the physical parameters of the gas are 
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approximately equal for all parts of the spherical region, the inte- 
gration can be reduced to a more simple expression: r v (L) w a„L. 
So, the optical depth in this case can be written as: 



A 3 (2J' + 1) 
8tt (2 J + 1) Vi 



L. 



(12) 



where A is the wavelength, J is the rotational quantum number, 
xud is the relative abundance of the HD molecule, Ay j is the 
Einstein coefficient, and nj and nj are the baryon and J-th ro- 
tational level population numerical densities at epoch 2, respec- 
tively. As to nj (2), it is well known that at the low density limit 
(tih < 10 3 cm -3 ), the population differs from the Boltzmann one. 
So, in the general case, when the gas is out of thermodynamic al 
equilibrium, the correct values of the ro-vibrational levels popula- 
tion should be calculated with the detailed balance equation: 



n vJ 53(W&_ MI /j/ + W°,^v'j') = 

v'Ji 

v'J> 



(13) 



where n v j is the population of the ro-vibrational level vj, 
W^j^ v /j/ and W£j_> v > j, are the probabilities of the radiative and 
collisional transitions, respectively. 

By using the expressions 1 1 21 and 1131 . we calculate the HD 
line optical depths corresponding to protohaloes of several masses 
at their turnaround redshifts (Table 1). The results corresponding to 
three ground rotational line transitions (1-0), (2-1) and (3-2) for 
the 3cr and 6<j protohaloes are plotted in Figs. 2 and 3, respectively. 
As one can see, the values of the optical depth for the HD molecule 
lines are rather large as compared with those reported for the LiH 
molecule (Bougleux & Galli 1997). Even the transition J' — J = 
3 — 2 reaches values as high as 10~ 9 — 10 -8 for the redshifts 
2 « 20 - 40. 

The estimates presented above correspond to the linear evolu- 
tion phase of fluctuations. The calculation of chemical kinematics 
evolution of HD molecule during the gas cooling and collapse in- 
side the virializing dark haloes is very difficult since it is linked to 
the thermal and hydrodynamical evolution. Under several assump- 
tions and simplifications, Galli & Palla (2002) presented results of 
HD abundances after the gas cooling inside virialized CDM haloes 
constructed according to Tegmark et al. (1997). The result of Galli 
& Palla is that the HD relative fraction, xhd, increases dramati- 
cally, approximately by a factor of 200, with respect to the initial 
one at the collapse epoch of the halo. In the following, we will use 
this result to get a rough estimate of the HD relative fraction in 
collapsed ACDM haloes and the corresponding t v . 

We have calculated above the evolution of the HD relative 
fraction, xhd, in the homogeneous pregalactic medium. For a 
given mass fluctuation va, we have obtained ihd at 2 ta . Let us as- 
sume that after the halo collapsed, at 2 co i, uhd increased only pro- 
portionally to the halo density increasing. Therefore, riHD(zcoi) = 
riHD(zta) x Pcoi/pta- Then, we assume that the gas cools rapidly 
and falls to the halo center in a dynamical time, 



t&y 



tt/2 



y/2GM/Rl 



(14) 



where i? co i is the halo radius. Now, according to the Galli & Palla 
result, we assume that at the epoch when the gas cooled and col- 
lapsed, 2 coo i, the HD relative fraction increased by a factor of 
~ 200, i.e. xhd(«cooi) ~ 200 ihd(zcoi)- The collapsing factor 
of the gas is assumed to be 10, so that the number density increases 
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Figure 2. Optical depths in the first, second and third rotational transitions 
of HD molecule in 3<r ACDM protohaloes of different masses reaching 
their maximum expansion at the redshifts shown in the abscissa. 




Figure 3. Same as in Fig. 2 but for 6<r ACDM protohaloes. 



by a factor of 1000. We further assume that the kinetic tempera- 
ture of the cooled gas is equal to the CMBR radiation (see Fig. 
10 in Galli & Palla). The populations corresponding to the main 
ground line transition (J' — J)=(l-0) actually depend weakly on 
Tfe. In Fig. 4 we show the HD (J' — J)=(l-0) line transition optical 
depths corresponding to 3 and 6<r haloes of several masses at their 
"cooling" redshifts. The radii used to calculate t u were l/10th of 
the corresponding halo radii R co i. For the redshifts of interest, i.e. 
before massive star formation (z > 20), the optical depths of the 
lines are significantly high, with values from several times 10 -2 to 
several times 10 , up to 2 ~ 50. 



3.3 Observational estimates 

As shown by Sunyaev & Zel'dovich (1970) (see also Dubrovich 
1977; Zel'dovich 1978; Maoli et al. 1994) the amplitude of sec- 
ondary anisotropies in the CMBR temperature due to Doppler ef- 
fect of resonance scattered photons in protoclouds with peculiar 
motions is: 
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Figure 4. Same as in Figs. 2 and 3 but for haloes already collapsed and 
at the corresponding redshifts when the baryons have supposedly cooled 
and collapsed inside them (see text). Only the case for the first rotational 
transition is showed. 



for t v « 1, where V p = V p (z) is the peculiar velocity of the 
protocloud with respect to the CMBR at epoch z, c is the light 
speed, and t„ is the optical depth at frequency v through the pro- 
tocloud at epoch z (to be the protocloud zt& in our analysis). The 
linear theory of gravitational instability shows that the peculiar ve- 
locity of every mass element grows with the expansion factor as 
V p (z) oc D(z)/(1 + z). An accurate approximation to this expres- 
sion for a flat universe with cosmological constant is (Lahav et al. 
1991; Carroll et al. 1992): 



V p (z) oc f(z 



9(z) E(z) 

5 (o) (i + z y- 



(16) 



where /(z) w Q.m{z) ' 6 , and the functions g(z), E(z) and SIm(z) 
were defined in §3.1. For the Einstein-de Sitter universe, eq. i 1 61 
reduces to the known expression V p (z) oc (1 + z)~ 2 . 

The peculiar velocity field at z = 0, traced by galaxies and 
galaxy clusters, has been measured in a large range of scales (for a 
review see e.g., Zaroubi 2002). For our problem, the velocity field 
traced by galaxy clusters is probably the interesting one because, as 
mentioned above, the first luminous objects in the universe formed 
in the densest, most rare regions, i.e. those that today are clusters 
of galaxies. Various data sets lead to different rms peculiar ve- 
locities for scales > 100h _1 Mpc. Here we will adopt a value of 
V p (0) ~ 650km/s, in agreement with results from Lauer & Pots- 
man (1994), Hudson et al. (1999,2004), and Willick (1999). This 
value is apparently a factor of ~ 2 larger than the one correspond- 
ing to ACDM numerical simulations for the scales of interest. 

By using eq. 1 1 61 normalized to V p (0) = 650km/s and the 
HD optical depths calculated in §3.2, we may estimate with eq. 
1151 the corresponding temperature fluctuation of the secondary 
anisotropics. In the range of z — 20 — 40, for the case of the 3<j 
protohaloes, we obtain that AT/T » 10 -10 - 4 10~ 12 in the 
first rotational transition, while for the 6a protohaloes, AT/T m 5 
10" 10 - 10" 10 .Inthe same redshift range, for the collapsed bary- 
onic clouds, whose HD abundances were estimated using results 
from Galli & Palla (2002), we calculate AT/T w 10 -6 - 10~ 8 
and AT/T w 2 10~ 5 - 10~ 6 for the 3 and 6cr cases, respectively. 

Now let us estimate the integration time required for the de- 
tection of the secondary anisotropics with millimeter and sub- 
millimeter facilities under construction as for example: the 50m 



LMT/GTM telescope in Sierra Negra, Mexico; the "Combined Ar- 
ray for Research in Millimeter-Wave Astronomy" (CARMA) in 
USA; the "Atacama Large Millimeter Array" (ALMA) in Chile. 
The observational (integration) time At can be estimated from the 
equation 



AT 



VAuAt 



(17) 



where Av is the bandwidth, T n is the detector noise temperature 
and AT is the amplitude of the temperature fluctuation calculated 
above. For the three facilities mentioned, T n w 40 — 80 K, and 
we may assume Av « 1. In the case of GTM/LMT telescope, its 
angular resolution limits the observability of several cases (see be- 
low). By using eq. <17> . for the 3<r and 6a protohaloes, we find that 
the integration times required to get an observable signal of the first 
rotational transition of HD molecules are too large to be observed. 
In the case of collapsed gas clouds (previous to star formation trig- 
gering) inside virialized 6a haloes, we estimate At < 10 s = 28 
hours for z < 40 and T n « 50 K. This integration time can be 
attained at different observational sessions of several hours each. 

Let us discuss another important parameter of the potentially 
observed protoclouds in molecular resonant line emission: their ap- 
parent angular size or diameter. Apparent angular size for an object 
of the linear size L is given by the Hoyle's formula: 



AO = 



L 



d A (z) 



(18) 



where cLa{z) is the angular-diameter distance. 

The angular-diameter distance for a flat universe with cosmo- 
logical constant is given by: 



<Ia{z) 



dz' 



(l + z)J H{z'Y 
where the Hubble parameter in this case is: 



(19) 



H(z) = H ov /n A f,o(l + z) 3 + «a, (20) 

and the radiation density term was neglected. Now eq. 1 1 8i can be 
written as follows: 



A0 - 



H L{l + z) 



c J (fi M ,o(i + z'Y + n A r 1/2 dz' 



(21) 



We use eq. 12 1 i to calculate the angular size A6 of 3 and 
6a protoclouds reaching their maximum expansion in the redshift 
range of 20 < z < 40 (corresponding roughly to mass ranges 
of 10 5 < M/Mq < 10 9 and 10 9 < M/M < 10 11 for the 
3cr and 6a cases, respectively). We calculate also the correspond- 
ing redshifted frequencies, vo, for the same three lowest rotational 
transitions of HD molecule considered in Figs. 2 and 3. Our re- 
sults are shown in Fig. 5, along with the expected coverage in this 
vq - A6 plane for the LMT/GTM and CARMA facilities (boxes). 
The ALMA facility covers the whole plotted plane. We also in- 
clude in this plot the curves corresponding to the cold clouds within 
collapsed 3 and 6a haloes for only the first rotational transition 
line (thin lines with square dots). These clouds could be resolved 
by ALMA and only partially by CARMA. It should be noted that 
high-peak halos are strongly clustered (Gao et al. 2005; Reed et al. 
2005). It may thus be more plausible to find a cluster of such ob- 
jects in a patch of the sky even if individual objects are very small. 
Moreover, the clustering increases the amplitude of the spectra- 
spatial fluctuations in the CMBR temperature. Figure 5 is rather 
general and can be applied to any other resonant lines of similar 
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Figure 5. Relation between the angular size AO of the 3a and 6<r pro- 
toclouds reaching their maximum expansion at 20 < z < 40 and the red- 
shifted frequency for three lowest rotational transitions of the HD molecule. 
Thin lines with squared dots are for 3<r (upper curve) and 6cr (lower curve) 
collapsed cold clouds in the same redshift range and only for the lowest tran- 
sition. The boxes indicate qualitatively the regions in which the LMT/GTM 
and CARMA will operate. ALMA will cover almost the whole plotted 
plane. 



frequencies in which the gas in protohaloes or haloes could be de- 
tected. 



4 DISCUSSION AND CONCLUSIONS 

The hope to detect signals from the end of the so-called 'dark age' 
has prompted a flurry of theoretical and instrumental activity. The 
observation of high-redshift protostructures before they collapse 
and star formation triggers will certainly help us to constrain struc- 
ture formation scenarios. One possibility of such observations is 
based on the spectral-spatial fluctuations in the CMBR tempera- 
ture produced by elastic resonant scattering of CMBR photons on 
molecules located in protostructures moving with peculiar veloc- 
ity. This kind of secondary anisotropies search for LiH molecule 
has been fruitless (de Bernardis et al. 1993). As Bougleux & Galli 
(1997) showed, the abundances of LiH molecule and their optical 
depths in the rotational lines are too small to produce detectable 
CMBR temperature fluctuations. 

In this paper we have investigated the spectral-spatial fluctu- 
ations in the the CMBR temperature due to elastic resonant scat- 
tering by HD molecules. For this molecule, the main contribution 
to the optical depth comes from the ground rotational transition 
(J' — J)=(l-0), because of small populations with high rotational 
levels of the HD molecule at early epochs. In the case of the LiH 
molecule, most of the contribution to the optical depth comes from 
rotational transitions between high J levels. The wavelenght at rest 
of the HD (J' — J)=(l-0) transition is 112.1 fim. Therefore, for 
z = 20 — 40 it should be observed at the wavelengths « 2— 4.5mm. 
For higher rotational transitions, the wavelengths are smaller. 

We have carried out calculations in order to estimate the HD 
optical depths and the corresponding amplitudes of the spectral- 
spatial fluctuations in the CMBR temperature produced by mov- 
ing protostructures. First we calculated detailed chemical kinematic 
evolution for HD molecule in the expanding and adiabatically cool- 



ing homogeneous medium, assuming standard BBN yields. We fur- 
ther calculated the optical depths in HD pure rotational lines for 
mass perturbations (protostructures) at their maximum expansion 
under the assumption that the HD relative fractions, xhd, within 
these protostructures are similar to the fractions calculated for the 
homogeneous medium at a given z. The spherical top-hat approach 
and the concordance ACDM cosmology were used to calculate the 
corresponding turn-around redshifts, sizes, and peculiar velocities 
for different masses and peaks. Once the optical depths for the first 
three ground rotational transitions were calculated, the correspond- 
ing fluctuations AT/T were estimated. We also estimated r„ and 
AT/T for cold gas clouds inside collapsed CDM haloes, but under 
several assumptions and using results on the HD relative fraction 
from a crude calculation by Galli & Palla (2002). 

The main results from our study are the following: 

• The relative fraction of HD molecule, ihd, increases in the 
homogeneous expanding medium by a factor of w 400 since z — 
300 to low redshifts (Fig. 1). At redshifts 40 - 20, xhu ~ 4 - 
5 10 -10 . These fractions are expected to strongly increase during 
the cooling and collapse of the gas within CDM haloes. 

• The optical depth of the HD first pure rotational line in high- 
peak protostructures at their maximum expansion increases with 
time. For the second and third lines, the corresponding t„'s at- 
tain a maximum at high redshifts and then strongly decrease for 
lower redshifts. The ranges of r„ values of the (1-0) transition 
lines in 3 and 6a protostructures at their maximum expansion in the 
20 < 2 < 40 redshift interval are approximately 10~ 7 — 6 10~ 9 
and 8 10~ 7 — 2 10~ 7 , respectively (Figs. 2 and 3). For the crude 
model of cold gas clouds inside collapsed 3 and 6a CDM haloes, 



the ranges of r„ in the same redshift interval are 8 10 



10" 



and 2 10 -2 — 1.3 10~ 3 , respectively (Fig. 4). The optical depths of 
HD rotational lines in protostructures are much higher than those 
of the LiH molecule. 

• The ranges of redshifted frequencies and typical angu- 
lar sizes in the sky of spectral-spatial fluctuations due to HD 
molecule rotational lines in high-peak protostructures and gas 
clouds in collapsed haloes fall partially within the observational 
windows of submillimeter telescope facilities under construction, 
as GTM/LMT, CARMA and ALMA (Fig. 5). The range of red- 
shifts under consideration was 20 < z < 40; after this epoch the 
first baryonic structures are probably already in site with strong 
emitting stellar sources. The observational search of the spectral- 
spatial fluctuations studied here will be important for testing mod- 
els of cosmic structure formation as well as the primordial abun- 
dance of deuterium predicted in standard and non-standard BBN 
theories. 

• For the ACDM scenario, the amplitudes of the spectral- 
spatial fluctuations produced by HD molecule at 20 < z < 40 
in protostructures, even as rare as 3 — 6a density-peaks, are too 
faint for reaching the flux sensitivity of submillimeter/millimeter 
facilities under construction. For the case of cold gas clouds within 
collapsed ACDM haloes emerging from ~ 6a density-peaks, the 
estimated amplitudes, in particular for the lowest rotational tran- 
sition line, are much larger and could be detected, but their esti- 
mated angular sizes in the sky are too small for telescopes such as 
GTM/LMT, but possibly clusters of such objects could be resolved 
(high-peak halos are strongly clustered). The spectral-spatial fluc- 
tuations predicted for very rare high-redshift cold clouds in the 
ACDM scenario will fall within the observational capabilities of 
other facilities as CARMA and ALMA. For non-standard BBN 
theories, the [D/H] abundance can be much larger than the used 
here. For example, in inhomegeneous BBN models, [D/H] is esti- 
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mated to be ten times larger than in the standard BBN (e.g., Lara 
2005). Therefore the HD line optical depths presented here could be 
a factor of ten higher, making already detectable the corresponding 
spectral-spatial fluctuations produced at the maximum expansion 
of protohaloes. 

A more detailed study of spectral-spatial fluctuations due to 
HD molecules in collapsing gas clouds as well as the inclusion of 
alternative BBN models is necessary in order to get more precise re- 
sults and predictions than the ones presented here. In this paper we 
have carried out preliminary calculations that show the viability of 
using HD molecule lines to search for signatures of early structure 
formation in the universe and/or to constrain cosmological theories. 
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